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Abstract 

Photon distribution amplitudes up to twist four are calculated within the nonlocal chiral quark 
model with a simple pole ansatz for momentum dependence of the constituent quark mass. Cal- 
culations are performed using modified electromagnetic vector current in order to satisfy Ward 
identities. Quark condensate and magnetic susceptibility of the QCD vacuum entering definitions 
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results are analytical up to the numerical solution of certain algebraic equation. 
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I. INTRODUCTION 



In the present paper we calculate photon distribution amplitudes (DA) within a low 
energy nonlocal model based on the instanton model of the QCD vacuum. There are seven 
different photon DAs corresponding to the Dirac structure probing the photon and to the 
light-cone twist (here we follow closely definitions of ref. However in reality only one 



of them, twist 2 tensor DA, can be accessed experimentally in hard exclusive processes 

mm 

(for experimental overview see ref. Q]). Higher twist amplitudes are suppressed 
in hard processes and vector twist 2 DA decouples for real photons. However, the interest 
in the remaining photon DAs is not purely academical. They are normalized through low 
energy observables such as quark condensate - {i'i'), magnetic susceptibility - Xm and 
mixed quark-gluon condensate - /s^ that are of importance for our understanding of the 
properties of the QCD vacuum. Only the calculation of all photon DAs provides a test of 
the whole approach and may prove its consistency. 

In the present work we employ a nonlocal generalization of the semibosonized Nambu 
Jona-Lasinio (NJL) model with momentum dependent constituent quark mass M{p) (which 
will be denoted as Mp) which follows from the instanton model of the QCD vacuuin [7 
This model in the present version has been previously used to calculate pion y, Um, 

nnn n ^ . 

and kaon [UJ, [IJ, ll3| distribution amplitudes, two pion DAs |15l|. p ion-to-photon transition 
distribution amplitudes jl^ and also twist 2 tensor photon DA jl7| (see also |^). However a 
complete analysis of all seven photon DAs has been, to the best of our knowledge, conducted 
only in ref. jl^ in a model similar to ours with, however, results that in some respects are 
different than the ones obtained in the present paper. 

One of the obvious problems that arises when one considers momentum dependent fermion 
mass is the nonconservation of the naive vector current containing only Dirac matrix. 
There are many proposals in the literature how to extend vector current to satisfy elec- 
tromagnetic Ward identities. None of them is unique, since current conservation does not 
fix the transverse part of the modified vertex. One of the simplest extensions of this type 



20|, 121| and employed also in ref. [l^, consists in 



discussed already many years ago in refs. 
the following substitution 

Extension ([T]) follows from the assumption concerning both Ward identities and analytical 



structure of the modified vertex that is required to match perturbative expansion. In prin- 
ciple one could add to ([T]) terms proportional to r'^ where r - p = 0, and the Ward identities 



would be satisfied. In refs. 



22 



23j and also in [24\ another modification has been advocated; 



here the ambiguity is connected to a freedom of choosing the integration path that defines 
the nonlocal vertex. In view of this ambiguous situation we have decided to use the simplest 
possible extension of ref. {2^ given by eq. ([I]). 

Unlike the pion or the p meson the photon has a dual nature being both point-like and 
composite at the same time. Therefore in order to calculate photon DAs that describe 
nonperturbative quark-antiquark structures, one has to subtract the perturbative part. In 
order to avoid ambiguities this procedure has to be well defined. In our case, since we work 
in the chiral limit, we subtract the perturbative part only from these photon DAs that are 
nonzero for massless quarks. These are: vector and axial DAs which are also UV divergent. 
Therefore the subtraction of the perturbative part renormalizes these DAs ensuring their 
finiteness. At the same time it introduces a term proportional to ln(— P^) that develops 
imaginary part for positive virtualities. This is the reflection of the fact that the photon 
can decay into free massless quarks in the chirally even channels. Throughout this paper we 
plot DAs both for negative and positive photon virtualities, in the latter case we take just 
the real part if the imaginary part exists. We also display momentum dependence of the 
pertinent decay constants that are characterized by dimensionless functions Fxy^AiP^) for 
brevity referred to as form factors. 

Finally let us make a technical remark concerning loop integrals over d'^k with = uP^ 
component fixed by the 5 function. Such integrals, depending on the tensor structure may 
contain "singular" pieces proportional to 5{u) and 5{u — 1) or even the derivatives of the 5 
functions. This is the case for higher twist photon DAs only. We discuss this in more detail 
in sect. HVl here we just want to point out that higher twist DAs are in fact distributions 
rather than ordinary functions. Quite importantly, the 5 function contribution is always 
accompanied by a regular piece that together with the 5 piece integrates to zero for any P^. 
This allows to define a properly normalized regular part and a singular part of DA of zero 
norm. 

In the next section we introduce kinematical variables and define photon distribution am- 
plitudes. In sect.milwe describe the main features of our model specifying the ansatz for the 
momentum dependence of the constituent quark mass. We fix model parameters requiring 
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that the experimental value of the pion decay constant = 93 MeV is reproduced. To this 
end we use model formula given in eq. (|T8|1 . Next, in sect. HVl we describe techniques used to 
calculate loop integrals with momentum dependent constituent quark mass. We pay special 
attention to Lorentz invariance and show how the end point singularities proportional to 
the Dirac 6 functions arise. Main results are presented in sect. |Vl First we calculate dimen- 
sional constants entering definitions of the DAs (I7D-([9]), namely quark condensate, magnetic 
susceptibility and yet another constant, called f^-y. We obtain numerical results that agree 
with the "experimental" values known from phenomenology. Finally in sections IV Al and 
IV Bl we present our main results plotting different DAs and discussing their properties. 

Our results can be briefly summarized as follows. Leading twist amplitudes are rather 
insensitive to model parameters, whereas higher twist amplitudes exhibit much stronger 
dependence, moreover some of them contain 6 functions. We also show the influence of the 
nonlocal part of the photon vertex ([ID on the shape of photon DAs. For some DAs it is rather 
unimportant, for the other ones it is absolutely crucial. More discussion and comparison 
with other models is given in sect. IVTl Technical details are collected in appendices. 



II. DEFINITIONS AND KINEMATICS 



Photon distribution amplitudes are deflned through matrix elements of the nonlocal 
quark-antiquark billinears between vacuum and one photon state. Quark operators are 
assumed to be on the light cone separated by the distance 2A. In the following we use the 
light-cone coordinates deflned by two light like vectors ra^ and such that = (1, 0, 0, — 1) 
and = (1,0,0,1). In this basis any four-vector can be decomposed into and v~ 
components 

v'' = t;+Y + w"Y + <. (2) 

Scalar product can be written as 

We shall work in the system where the photon momentum is expressed as 

2 P+ 2 ^ ' 
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Decomposition of the polarization vector reads 

= V s~ where £^£±^ = — • Cl = — 1- (5) 

Since P • £ = we have the relation 

e- = ^£+. (6) 

For real photon we obviously have — e ■ n — ^ and consequently £~ = as well. 
Depending on different tensor nature of the bilocal operators, we can define tensor 

(0 I (An) (-An) | 7 (P, e)) = (V^V^) (P^) 

I [elvP - £^n") P^xm du e^«^^" 0^ («, P^) 

(^«^/3 _ ~/3^a^ ^+ / due'^^^^ iPt{u,P^) 

P Jo 

^ (£>'^ - sin-) du e^«^^" hT («, P^) ^ (7) 

vector 

(0 (An) (-An) \l{P,e))^ eh,Fy (P^) 

2 Jo 

+ £^[ I' due'^^''^iPv{u,P^) 
Jo 

--^^n^s+ C due'^^^'^hviu.P^) 
2{P+f Jo ' y 

and axial vector 



(0 1^ (An) 7^5^ (-An) | 7 (P, £)> = ^e/s^P^ (P^) 



io 



(9) 



distribution amplitudes. For compactness we used notation ^ = 2m— 1 where u is longitudinal 
fraction of the quark momentum and dropped Wilson lines [—An, An] that ensure gauge 
invariance of the nonlocal operators. In the light-cone gauge A-n — and hence [—An, An] = 
1. Moreover, since we use an effective model where gluonic fields are integrated out, Wilson 
lines corresponding to gluon fields never appear. 



fl and 0, 



Our definitions follow closely those of refs. [Ifl and [18|, however we need only one 
P^— dependent dimensionless form factor for each tensor structure: Fj'(P^), Fy (P^) and 
Fa (P^) where subscripts T, V and A stay for vector, tensor and axial vector, respectively. 
Constant Xm is the magnetic susceptibility of the quark condensate (V^V^), and /s^ corre- 
sponds to the axial mixed quark-gluon condensate. They provide natural mass scales for 
distribution amplitudes. Analytical expressions for {i'i'), Xm and /s^, and for the form 
factors can be retrieved from the matrix elements of local operators: 

(0 1^(0) a'^"^ (0)1 7 iP,e)) = le {^^) XmFr (P') {e'^P'' - e^P") , (10) 
(0 1^ (0) j^i; (0) I 7 (P, e)) = efs.Fy {P') £^ (11) 

= efs-yFA (P') e^^^f^e'^P'^n^ (12) 



^(0|^(-An)7^75^(An)|7(P,£)> 



A=0 

Equations ((ZD"© define photon distribution amplitudes that can be classified according 
to the kinematical light-cone twist. We have distributions of twist-2: 0t, (pv, of twist-3: 
^T, ^y, i^A and of twist-4: /i-r, hy. This can be easily seen by inspecting eqs. (I7I)-(l9D, since 
the twist counting actually reduces to counting the powers of P"*". Notice that in the case 
of axial distribution the power of P"*" equals 1, what would correspond to twist-2, however 
additionally there is a path stretch A with inverse mass dimensionality that makes ipA to 
have twist 3 rather than 2. 

Constants Xm, and /s^ are chosen in such a way that the following normalization 

conditions are satisfied 



/ <l)T{u,P^)du = l, [ ^T{u,P^)du = XmP\ f hT{u,P^)du = XmP' 

Jo Jo Jo 

I (pv {u, P^) du=l, I ipv {u, P^) du = l, I hv {u, P^) du = 1 

Jo Jo Jo 

[ ipA {u, P^) du = 1. 
Jo 



(13) 
(14) 



(15) 

Note that due to the conservation of vector current Fv{0) = 0. On the other hand 
Pt(0) = 1. Normalization of the axial form factor Fa{0) is arbitrary and depends on the 
dimensional constant used in definition (191). 
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III. NONLOCAL CHIRAL QUARK MODEL 



In order to calculate relevant matrix elements in the low energy domain we shall use 
effective action based on the instanton vacuum theory Its main feature is momentum 
dependent constituent quark mass 

M{k) = MF^{k) (16) 

appearing due to the chiral symmetry breaking. This dependence enters not only into 
propagators, but serves as a nonlocal quark-meson coupling as well. For zero momentum 
M (0) = M is of order of 350 MeV, while for /c ^ oo constituent quark mass vanishes 
M{k) 0. 

One has to remember that the semibosonized NJL model, although devised to describe 
chiral physics of Goldstone bosons, has been widely used to incorporate baryons as chiral 



solitons both in local (for review see e.g. ref. [25|) and non-local [26|| cases. Generally the 
results of these studies show that the soliton ceases to exist for too small constituent quark 
mass M. The critical value of M depends on the details of the given model, however it is of 
the order of 300 MeV or a bit less. Typical values of M that fit well the hyperon spectrum 



may be as high as 420 MeV [27] . In order to investigate dependence of photon DAs on M 
we use three distinct values of M: 300, 350 and 400 MeV. 

Due to the momentum dependence of the quark mass, the naive vector current j'^ = 
ip'jf^ip violates electromagnetic Ward identities. In order to fix this deficiency new nonlocal 
terms have to be added to j'^. As already discussed in the Introduction there are several 
ways of constructing extensions that make the current conserved. In the present paper we 
use the simplest possible generalization of the vector current, 13 replacing 7^ by an 
effective vertex of ([I]) . One can easily check that electromagnetic Ward identities are satisfied 
when ([T|) is used instead of 7^^. Although eq. ([I]) introduces an extra pole inside Feynman 
amplitudes, its residue, as we will explicitly show, is zero due to the mass difference in the 
numerator. This generalization of the vector current has been widely used in the literature 



also in the context of the photon DAs 



18|. 



Expression for the form factor Finst [k) within instanton vacuum model is known analyti- 



cally in Euclidean space and is highly nontrivial 



7|. Therefore, in order to perform analytical 
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calculations directly in Minkowski space, we use the following formula 



where A„ is cutoff parameter adjusted for each n in such a way, that the experimental value 
of the pion decay constant is reproduced. For transparency we shall skip subscript n and 
use A rather than A„ in the following. 

Equation (flTll reproduces reasonably well original shape Fjnst (k) when continued to Eu- 
clidean momentum. It should be however pointed out that expression (fTTl) does not follow 
the exponential asymptotics of Finst (k) [7|. Parameter n is introduced in order to check 
sensitivity of our results to the shape of F {k). 

In order to fix the model parameter A we use the following Euclidean expression for the 



weak pion decay constant 



22|: 



F^ = ^ r dk' e - fcjM (fcj) M' + fc^M- jklf 

where M'(fc|) = dM{k^)/dk^. Notice that this formula differs from the Pagels-Stokar 



formula of ref. 



n 

28l |. It has been also obtained in ref. [10(| from the PC AC relation 



in Minkowski space. Using experimental value F^, = 93 MeV and (1171) we obtain cutoff 
parameters listed in Table [T] for several choices of the constituent quark mass M and n. 
Analytical expression obtained within the present model is given in appendix [Bl We remark 
at this point that the cutoff parameter A should not be confused with a typical scale of the 
model, which for the instanton model is about 600 MeV. 



IV. LOOP INTEGRALS WITH MOMENTUM DEPENDENT MASS 

In this section we present a brief sketch of our calculations underlying the most important 
steps. Further technicalities are relegated to appendices. In order to calculate the DA of 
interest - denoted generically as f{u) - we have to invert formulae (I7I)-(l9|) by contracting 
them with appropriate 4- vectors and by performing Fourier transform in A. This results in 
the following formulae 

— ?> 4P+/V r d^k 
f{u) = ^^^^1^ j AA. Tr{k, k - P) 6{k^ - uP^) (19) 
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M = 300 MeV 


n = 1 


A = 1016 MeV 


n = 5 


A = 2385 MeV 


M = 350 MeV 


n = 1 


A = 836 MeV 


n = 5 


A = 1970 MeV 


M = 400 MeV 


n = 1 


A = 721 MeV 


n = 5 


A = 1704 MeV 



Table I: Numerical values of the model parameters obtained using Birse-Bowler formula for pion 
decay constant F-^. 

where C is the constant obtained by the contraction, T is the contracted tensor structure 
defining given amplitude and 

r = ixr jr- -e^Y{k,k- P)— \, -I (20) 

stands for the Dirac trace. Note that in the case of axial DA, because of A standing in the 
l.h.s. of (HD 

/(«) = \ + Mmu) - i^A5{u - 1)) . (21) 

Since some of the integrals can be UV divergent we shall work in D = 4 — 2e dimensions. 

Previous calculations using present nonlocal model were done by integration in the light- 
cone coordinates, with special care concerning the integration contour to ensure analyt- 
icity in A„. Here we present another method of performing such integrals based on the 
a-representation for the propagators. It is especially useful in the case of integrals appear- 
ing in higher twist distributions, because of the end point delta-type singularities, which are 
cumbersome to treat by integration in the light-cone coordinates. 

The above complication can be well illustrated by considering a loop integral of the type 
ffT9|) for the numerator involving k^. If not for the 6 function it would have been proportional 
to P^, but because fc"*" — uP^ = n- k — un-Pwe have 

— 4t / , -6(k+ - uP+) = A(u)P'' + B(u)n>' (22) 

{27if{k''-Ml+ie){{k-Pf-Ml_p + ze)^ ' I i 
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where M is some scalar function involving n, e and P. There is an obvious condition following 
from Lorentz invariance 



j duB{u) = 0. (23) 



However, as mentioned above and as shown explicitly in appendix O function B{u) contains 
both regular piece and the piece with delta functions: 6{u) and 6{u — 1). Only the sum of 
both contributions integrates over du to zero. Note that this cancelation occurs for any P^. 
Since n ■ n = and e± ■ n = 0, the delta functions contribute only to the integrals of the k~ 
component. The integrals with tensor structure k'^k" are even more complicated, since they 
involve derivatives of 6 functions. 

As it was already discussed in ref. momentum mass dependence given by (ITTl) intro- 
duces a set of poles, whose positions depend on parameter A. To this end it is convenient 
to introduce dimensionless scaled variables 

K = k/An, p = P/An, r = M/An (24) 

and to define 

zi = {k — pf — 1 + ze, Z2 = — 1 + ie. (25) 

Then the loop integral involving two propagators, like the one in. eq. (I22l) . is transformed 
into 

where the numerator jV^zi, Z2, n ■ n, n ■ n, k ■ e±) depends on the DA considered. Here 

4n+l 

G (z,) = zt+' + zt -r'=l[{z,- V,) (27) 

j=i 

corresponds to the propagator with momentum dependent mass (for n = it reduces to 
the ordinary propagator in scaled variables) where the complex numbers rij are roots of 
polynomial G to be obtained numerically. 

Next we decompose the inverse product of G (zi) G (^2) into a sum of simple poles 

^M^N 4n+l „M AT 

E f^f^ l..^ML-r.A M,N<An (28) 
'■J 

with 



^(^1)^(^2) ,j^/''izi-V^di^2-Vj) 

f =Il ^r—r- (29) 
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In this way integral (1261) is reduced to the sum of contributions involving two propagators 
only. It is convenient to use the a-representation (exponential Schwinger representation) for 
the product of propagators since also the 6 function in fl26l) can be written as an exponent. 
Further calculations are rather standard and are summarized in appendix O As a result 
we obtain analytical expressions given as sums over roots rji. Certain simplifications occur 
when we use the following identity (which is true for any set of complex numbers {rji} not 
only for the solutions of G{rii) = 0): 

.r I iV < An, 

E f^^^ = (30) 
i=i I 1 = An. 

The proof of fl30l) and other useful identities can be found in appendix [Al 

Some of the loop diagrams discussed above are UV divergent and require renormalization. 
This results in the subtraction of the perturbative part which is uninteresting from the point 
of view of the hadronic component of the photon. To illustrate this problem consider loop 
integral (l26l) with Af = 1. Performing d^K integration gives 

4ra+l 



4n+l . . 

E /^^^ /i^r \- - In [1 - + nr^, + nup^] j (31) 



167r2p+ V A2 ^ 

with u = u — 1. Because of (l30ll the coefficient of the 1/e pole is equal to 1. For M involving 
negative powers of rji (like N = for example) the coefficient of the pole is and no 
subtraction is needed. For constant mass (n = 0) ^(2;) = z + 1 — r"^ and for zero mass 
(current masses are zero in the chiral limit) there is only one solution of G{z) = 0, namely 
rji = —1. Hence the perturbative part of the loop integral (l26l) reads 



'^--i6;W(^)'(7-"'[''*"''0- '''' 

This result can be of course obtained by standard techniques for M = 0. Renormalization 
in the MS scheme proceeds by subtracting the pole only. Here we subtract full perturbative 
contribution and go back to D = 4 (e = 0) dimensions which gives 

4n+l 



167r2p+ 



l-urii + ur]j + uup'^] ^^^^ 



uu p2 
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where we have used again (l30ll . Note that subtraction occurs only for terms which do not 
involve or Mk_p, and these terms are always UV divergent. In other words in the chiral 
limit all perturbative photon DA are either UV divergent or identically zero. 



V. PHOTON DAS IN NONLOCAL MODEL 

Before we proceed with photon DAs and systematically present our results we have to 
fix numerical constants appearing in the definitions ©-(jQl). We have already introduced 
the expression for pion decay constant (|T8ll which is used to fix model parameters. Next 
we consider the quark condensate given as the trace of the quark propagator which reads in 
Euclidean metric 

which in our model turns out to be simply 

AT A/f A2 

(V^V^) = E /^^^ (1 + 1^ (1 + • (35) 

i=l 

Numerical values obtained from this formula coincide with those of ref. jsLif we use model 



parameters corresponding to obtained from the Pagels-Stokar formula 



28|. 



The formula for magnetic susceptibility Xm in the nonlocal model used in ref. 



(with M'{k^) = dM{kE)/dk^) reduces in our case to 

N M ^""^^ 

= E f^f^vt (1 + V,) [VT + 2n (1 + Vj) vf-'] (37) 

I (log (1 + V^) - log (1 + m)) + I 

where e^j is for i = j and 1 otherwise, while 5ij is Kronecker delta. Numerical values of 
(?/'?/') and Xm for the present set of model parameters are listed in Table HTl Note that in fact 
we do not have to use (1371) to calculate Xm since it can be retrieved from the normalization 
condition of 0t(m). Numerical values of Xm obtained both in ways agree proving consistency 
of our calculations and definitions 

To calculate /s^ we have used Euclidean formula from ref. [isl : 
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M = 300 MeV 


n = 1 


(^^) = -(277MeV)^ 


X„^ = 2.30GeV-2 


n = 5 


(^V') = -(230MeV)^ 


X„^ = 3.75GeV-2 


M = 350 MeV 


n = 1 


(^^) = -(253MeV)^ 


X„^ = 2.85GeV-2 


n = 5 


(^^) = -(208MeV)^ 


X„^ = 4.71GeV"2 


M = 400 MeV 


n = 1 


(^^) = -(236MeV)^ 


X„^ = 3.34GeV-2 


n = 5 


(^^) = -(192MeV)^ 


Xm = 5.58GeV-2 



Table II: Numerical values of the quark condensate {tpip) obtained using model parameters from 
Tab m and magnetic susceptibility Xm used in the calculations. 



M = 300 MeV 


n = 1 


= -0.0095 GeV^ 


n = 5 


/g^ = -0.0093 GeV^ 


M = 350 MeV 


n = 1 


= -0.0095 GeV^ 


n = 5 


= -0.0092 GeV^ 


M = 400 MeV 


n = 1 


= -0.0094 GeV^ 


n = 5 


= -0.0091 GeV^ 



Table III: Numerical values of /s^ obtained using model parameters from TabUl 



which in our model transforms into 

/37 = ^ + (39) 

i=l 

Numerical values of f^^ are listed in Table lllli 
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Phenomenological values of {ipip/, Xm and /s^ are well known only for the quark con- 
densate: approximately —(250 MeV)^ [29|| at low momentum scale. This value is still used 
in more recent phenomenological applications Magnetic susceptibility is still a subject 
of large phenomenological uncertainties. Different estimates are nicely summarized in ref. 
30| where it is shown that Xm — 2.5 -h 5.5 GeV~^ with some preference to the values around 
4.3 GeV~^. Finally the value of fs^ obtained in different low energy models, as discussed 



in ref. [l8|, is negative and of the order of —0.004 GeV^. Our values are here factor of 2 
smaller (~ 0.0094 GeV^), however they are almost insensitive to actual model parameters. 
On the other hand magnetic susceptibility is quite sensitive to M and n (see eqs. (fT6ll and 
fflTl) ) remaining, however, within the range of acceptable phenomenological values discussed 
in ref. Similarly {ijjtjj) varies with M and n, however, for the preferred value of the 

constituent quark mass M = 350 MeV it is quite close to the phenomenological estimates. 
From this point of view our model satisfactorily describes low energy observables relevant 
for photon DAs. 

A. Leading twist distributions 

1. Tensor photon DA 

Tensor twist-2 amplitude has been already discussed in refs. 0| and also jsl in a model 
with the local vertex only. Here we extend discussion to off-shell photons and calculate the 
correction appearing due to the modified vertex ([I]). In the case of leading twist tensor DA 
we obtain the following expression using local current: 

0(0) ( p2x ^4Ay^_^ f ^s{k-n-un-P) 
^ ' (^^> XmFr (P2) J (2vr)^ ^ ' 

u Ml. — u Mu_p , , 
— f40) 

(P-M| + ze)((A:-P)^-M|_p + ze)' ^ 

where u = u — 1. Notice, that special choice of the contour described in Q| allows for passing 
to Euclidean space. Therefore we can use Schwinger representation for scalar propagators 
and proceed in the spirit of [l^- The result reads: 

(41) 
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for < M < 1. 

Now let us consider the part coming from the nonlocal part of the vertex ([T]) . It is given 
by the integral 

(Mfc - Mfc.p) (ex ■ A;^)' 



(fc2 - M| + ?e) ((fc - P)' - M|_p + le) {2k ■ P - + zfi) 
Notice that there appears additional denominator 2k ■ P — P^ + ifi, where prescription 
introduced at this stage is completely arbitrary. However, as already explained in sect. [111 
the residue of this pole is zero so that it does not contribute to the amplitude irrespectively 
of the sign of /i. Therefore in the following we shall always omit contribution of this spurious 
pole. After performing the integration as described in appendix ICl we finally obtain 



(l + uup^ — urji + M?7j) In (l + uup^ — urji + urjj) 

for < M < 1. The full twist-2 tensor photon DA is given by 

A _ ^(0) , ^(1) 



+ <P't'- (44) 



We plot this function in fig. [H for several values of photon virtuality and constituent quark 
mass. Notice that the nonlocal part of the quark-photon vertex is small and the full ampli- 
tude is almost equal to the local one. The resulting DA is almost fiat for real photons and 
does not vanish at the end points. 

Tensor form factor is shown in fig. IV A 1[ It can be in principle calculated by analytical 
integration, which has to be performed carefully because of the complex numbers under 
logarithms. 

2. Vector photon DA 

Calculation of vector twist-2 amplitude proceeds in a similar way. After performing the 
traces we get 

-iAP+ N, f d^k + ^) Hk-n-un-P) 



^^^""^ h,Fy{P^)e+l {27^)D - Ml + ie){{k - Pf - Ml_p + le)- ^^^^ 
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b) 

M = 350 Me V, n = 1 <PtM 



.- =-03 GeV^ 



0.0 0.2 0.4 0.6 0.8 1.0 



1.0 
0.91- 



= 0, « = 1 




M =400 MeV 

M =350 MeV 

M =300 MeV 



0.0 0.2 0.4 0.6 0.8 1.0 



(f>T(u) 



M = 350 MeV, = 




M = 350 MeV 
n=l, P^=0 



J- u 




Figure 1: Leading twist tensor photon DA for: a) different photon virtuahties and fixed M = 
350 MeV and n = 1, b) different M and fixed n = 1 and = 0, c) different n and fixed 
M = 350 MeV and = 0, d) decomposition into contributions corresponding to local (dashed) 
and non-local (dotted) parts of the vector vertex for M = 350 MeV, n = 1 and P^ = 0. 



a) 

Fr{P-] 



n = 1 



M =300 MeV 

M =350 MeV 

M =400 MeV 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 



-p2[GeV2] 



M = 350 MeV 
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Figure 2: Tensor form factor for: a) fixed n = 1 and different M, b) fixed M = 350 MeV and two 
different n = 1,5. 
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where Ty^ and Ty^ stand for traces corresponding to local and nonlocal parts of the photon 
vertex respectively: 

ri°^ = e+{MkMk-P + kl- P^u) - {ei_ ■ kA.)P^{u + u), 

^ (46) 



^(1) _ (A4 - Mu-p) {uMu + uM,^p) 



£^\k- - - 2 (£l • kA_) 



Note that single powers of (el ■ /c^) integrate to zero. 

In the case of twist 2 vector DA we have to subtract the perturbative piece corresponding 
to M {k) — 0. Then, the contribution to the vector photon DA coming from the local part 
of the vertex consists of two parts 

= 4.^/3,;^(P2) E /^/.- i^'^t'^r (-^^ - - 1) + ^'^r^f ) 

In (l + uup^ — urji + ur]j) . (47) 

and 



The addition coming from the nonlocal part of the current can be conveniently split into a 
sum of two contributions 

In ^1 + uu]?' — ufji + ufjj^ , (49) 



In ^1 + uup^ — UTji + urjj^ . (50) 



Notice that subtraction concerns only 4>y'^^ part since (py^ is always proportional to mass. 

In the tensor case the only effect due to the local vertex is a small change in the shape of 
the distribution. The situation is different for vector DA. When we use local current only, 
the vector distribution alternates in sign (recall that leading twist DAs have probabilistic 



interpretation). Only when we include the nonlocal part of the vertex, contributions 0^°'"^ 
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Figure 3: Leading twist vector photon DA for: a) different photon virtualities and fixed M = 
350 MeV and n = 1, b) different M and fixed n = 1 and = —0.3 GeV^, c) different n and fixed 
M = 350 MeV and = —0.3 GeV^, d) decomposition into different contributions, as described in 



the main text; notice the exact cancelation of < 
For positive P^ we show only real part of DA. 



,(0,a) 



and 



following from the gauge invariance. 



and 



V 



cancel exactly for any P^. This is explicitly shown in figlUd. As a consequence 
(t)v is effectively the sum of 0^ '^^ and . We plot this function in fig. [3] for different sets 
of model parameters. 

Furthermore, since both 0y'*'' and 0^ ''^ are explicitly proportional to P^, normalization 
conditions ffT4|) require that -FV(O) = as it should be in accordance with the conservation 
of the vector current. This condition would be violated if not for the nonlocal part of the 
photon vertex. 
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a) b) 




Figure 4: Vector form factor for: a) fixed n = 1 and different M, b) fixed M = 350 MeV and two 
different choices of n = 1, 5. Notice that the form factor vanishes for zero virtuahty as required by 
vector current conservation. 

B. Higher twist distributions 

For higher twist distributions we encounter an additional difficulty. As shown in appendix 
ICl it turns out that they are in fact generalized functions. Due to k~ = k ■ h occurring in 
the numerator {k~^ = k ■ n is fixed - see delta function in (liOl) ). additional end point delta 
functions appear. They are crucial for Lorentz invariance of the integrals and in consequence 
for the correct normalization of the distributions. These singularities were already discussed 
in ref. LJ- 



1. Tensor DAs 



We start with tensor DAs. For twist 3 tensor amplitude we have 

4iV, P+ r d^k (^f^ + ^{k-n- uP+) 



(51) 



Ft{P^) e+ J (27r)4 (P - Mf + ie){{k - Py - M^_p + te) 

where T^°^ and T^^ stand for traces corresponding to local and nonlocal parts of the photon 
vertex respectively: 



^(0) 



— e 
2 



P+{Mk-Mk_p)[k- +u 



EL 

p+ 



Tr 



(1) _ I P+ (Mfc - Mk-p) 



k- + Uj^ 



p2 



^ \k - u— I - 2A;_L ■ e± 



(52) 
(53) 
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Note that even powers of k± ■ e± integrate to zero under d'^k±. 

Luckily in the case of ipriu) contributions involving k~ cancel out in the sum of (l52|l and 
fl53ll and the final result is the sum of two pieces 

NM 



(54) 



4n+l 



and 



P2 



„2n _ 2n 



In ( 1 + uup — urji + urjj j . 



(55) 



Note that ipr is proportional to the same normalization constant as (j)T times {xmP"^) which 
means that it decouples for real photons. 

In the case of twist 4 tensor amplitude the 5 function contributions do not cancel out. 
Performing Dirac traces we have 



hriu) 



(56) 



[i)^) FriP^)] (27r)4 {k^ - + ie){{k - Pf - M|_p + le) 

where again the contributions of local and nonlocal parts of the vector current have been 
singled out: 



R 



(0) 



e+{Mk-Mk-p) (e±-kA~i 



R 



(1) 



—I 



p+2 

P^ Mk - Mk-P 
P+2 k-+u^ 



{Mk - Mk-p)k- 



M, 



P+ 



e± ■ k_i 



1 - 2 ( £l ■ A;x 



(57) 
(58) 



Note that significant simplifications occur since single power of k± ■ e integrates to zero. 
Following the steps described in appendix [C] we finally arrive at the final formula for hpi^u). 
It is convenient to split it into 4 different pieces - regular local vertex contribution: 

4n+l 



NM 



In ^1 + uup^ 



UTji + ur]j 



(59) 



20 



second local vertex contribution: 

AT M 4n+l 

In ^1 + uup'^ — urji + urjj^ (60) 
that integrates to zero with 5-function contribution 

[I + r;,) In (1 + r/,) 5{u - 1) - (1 + r^,) In (1 + r/,) (61) 

and hence does not contribute to the normalization, and the contribution corresponding to 
the nonlocal part of the photon vertex: 

N M 471+1 

^ ^ In ( 1 + uup'^ — -ur/j + wij ] . (62) 

The delta contribution /i^''^'^'*^'* can be rewritten using the expression (l39ll for /s^ and the 
identities given in Appendix [Al 

= {Kn - 1) + (63) 

Note that the fact that the sum of (l60ll and (IBB integrates over du to zero is a consequence 



of Lorentz invariance discussed in sect. [IV] (see eq. fl23l) ). Therefore only /i^'"'' and /i^"*! 
contribute to the normalization condition (fT3l) given by (x^P^). If not for the 5-term /i^''''* 
would also contribute to the norm spoiling the normalization condition. 

Full results (with nonlocal current) for twist 3 V't and twist 4 Ut tensor distributions are 
show in figs. [5] and [6] respectively. 

2. Vector DAs 

In the case of higher twist vector DAs, and hy ^ calculations are basically the same, 
with the restriction that we have to perform subtractions similarly to the twist 2 case. For 
twist 3 amplitude we obtain: 

z4Ar,P+ /- d^k (^f + ^f) Hk-n-un-P ) 
h,Fy{P^)] {2nr {k' - Ml){{k - Pf - Ml_p) 



u 
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Figure 5: Tensor twist 3 photon DA for: a) different photon virtuahties and fixed M = 350 MeV and 
n = 1 (for = it is identically zero), b) different n and for fixed M = 350 MeV, c) different 
M and fixed n = 1 and = —0.3 GeV^, d) decomposition into contributions corresponding to 
local (dashed) and non-local (dotted) parts of the vector vertex. 



with 



(0) 



2(£"1 ■ k, 



k 



e^(e± ■ k±) ( k 



u- 



'p+ 



^(1 - 2u)P+k- - ^uP^ - MkMk-p. 



(65) 



In fact after integrating over the transverse angle the terms in the first line vanish. Next 



T, 



(1) _ (^"1 ■ k^) {Mi - M|_p) 



k- + Uj^ 



e^[k-- ) - 2{k± ■ e±) 



(66) 
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Figure 6: Tensor twist 4 photon DA (without end point delta functions) for: a) different photon 
virtualities and fixed M = 350 MeV and n = 1, b) various M and fixed n = 1 and = 0, c) 
two choices of n = 1,5 and fixed M = 350 MeV and = 0, d) decomposition into contributions 
corresponding to local (dashed) and non-local (dotted) parts of the vector vertex. 



Again only quadratic term {e±_ ■ k±Y survives integration over the transverse angle. The 
result can be split into a regular part coming from the local part of the vertex: 



(O.reg) 



u 



((1 + 2uu) rii-r,fP' + 2M' + (1 - 2u){^. - Vj)v-''vT^') , 



the part with delta functions also coming from the local part of the vertex: 

4n+l 



(0, delta) 
V ' 



U 



8vr^/3^Fy(P^) ^ L 



(67) 



(68) 
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and the nonlocal part: 

— ^ In (l + uup^ — UT]i + urij) (69) 

?7i -r]j 

The part with delta functions can be rewritten as 

^?''^"^^(«) = ['^(^ - 1) + ^(«)] , (70) 

where we used (l39l) and the identity A^r7f"(l + ?7i) = following from equation fl271) for 
zeros of G{z). 

Next we calculate twist 4 vector distribution amplitude hviu): 

iAN,P+ (P+) f d^k [^f +R^y^') 6{k-n-un-P) 

J 

where the traces read 

P+ (P+) 

4'^ = 1 ^ [Ml - Ml_p) - Mk (Mfc - Mk-p)] (73) 

The only additional complication is due to the second derivative of delta function - the 
details can be found appendix ICl eqs. (IC20l) - (|C22p . Results for the local part read: 



^^^""^ h,FyiP^) 2P2 J (2vr)4 (P - Mi)iik - Pf - M^p) ^^^^ 



rP = (k-y - ^k~ - - MkMk-P (72) 



(^2uu ri'^rifP^ - 2M' + r^f^r^f (2r/, - Vj + I - Mv^ - Vj)) A' (74) 
+ ^ p2 I (75) 



and 



delta) ( ) ^ ^ ^ ^1 _ 2M2 y /. In (1 + 



[(1 + r/i) [5(n)+5(^i-l)]+p25(«)] 
+ A' 2^ /. /i^^fr/f - In ^ J (1 + uup^ - ur^, + «r^,) 



+ (76) 
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and for the contribution coming from the nonlocal part of the photon vertex: 
'''^ ^""^ " 8n^h^Fv{P^) '^^P^ £ ^ ''''^^ ~ ''''' ^ 

2n 2n 

^^^4^ (r^. - V, - {2u - l)p') ( (r^r + r^f ) - / r^f ) , (77) 

= , '^,p,, E (1 + 1^ (1 + [^(«) + - 1)] • (78) 

OTT j-3-yl'v[P ) P ^ 

We note that 

thus it cancels with jj^ /ly' '^'''*'^^(u) as can be easily seen. 

Our results are shown in figs. [7] and [8] for twist 3 and twist 4 respectively. The magnitude 
of both distributions is growing unlimitedly when photon becomes softer (obviously distri- 
butions multiplied by the vector form factor remain finite). Notice however that for the real 
photon hv decouples. 

3. Axial DA 

We have only one distribution in the axial vector channel which is of twist 3. When 
inverting the definition ([9]), due to presence of A on the right hand side, we obtain the 
expression for the derivative of DA rather then for DA itself 



(w) = V^A («) + i'A (0) 5 {u) - i)A (1) 5 in) (80) 

_ . SAT, /■ dPk tJ^ , _ . 

' Fa{P^) J {2^f {k- - Ml){{k - Pf - Ml_,)' ^ - " ^ 

with 

Tf = (k~ - j^uj - e+P' [e^ . k^). (81) 

In the case of axial DA the nonlocal part of the vertex does not give any contribution. This 
is simply because the Dirac trace is equal to zero. To obtain ipA iu) one has to integrate 
fl80l) over du' from to m: 

u 

ijA{u) = jij'A{u')du'. (82) 
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Figure 7: Vector twist 3 photon DA for: a) different photon virtuahties and fixed M = 350 MeV 
and n = I, h) various M and fixed n = I and = —0.05 GeV^, c) two choices of n = 1, 5 and 
fixed M = 350 MeV and = —0.05 GeV^, d) decomposition into contributions corresponding to 
local (dashed) and non-local (dotted) parts of the vector vertex. 



Notice that the end point contributions cancel out and one might get an impression that 
tpA (u) is determined up to a constant. Fortunately we have at our disposal an independent 
formula for Fa{P'^) given by eq. (fT2ll (see also eq. (ID8I) in appendix[D|) and the normalization 
condition (|T5ll that fix the value of ^^(0) = ^a(I) 7^ at nonzero value. 

As in the case of vector twist 2 DA iPa{u) is UV divergent and requires subtraction of 
the perturbative part. The result splits into a regular part 

# '(^) = ;.o.^ y. h h Vtvr {{m-Vj)^' + {l- 2u)P') (83) 



1 + uup"^ — urji + urjj 



In 



uup^ 
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Figure 8: Vector twist 4 photon DA for: a) different photon virtuahties and fixed M = 350 MeV 
and n = f , b) various M and fixed n = 1 and = —0.05 GeV^„ c) two choices of n = 1,5 and 
fixed M = 350 MeV and = —0.05 GeV^, d) decomposition into contributions corresponding to 
local (dashed) and non-local (dotted) parts of the vector vertex. 

and the piece involving 6 functions: 

4?i+l 



i=l 



Note that rjf"' (1 + rji) = r^, and in virtue of ([39] 



(84) 



S{u) - 5{u - 1) 



(85) 



The form factor Fa{P'^) and ipAiu) itself is shown in fig. [9l We obtain the following values 
for axial form factor at zero momenta: Fa (0) ~ 0.77 for M = 350 MeV, Fa (0) ~ 0.79 for 
M = 300 MeV and Fa (0) ^ 0.75 for M = 400 MeV. 
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Figure 9: Axial photon DA for: a) different photon virtuahties and fixed M = 350 MeV and n = 1, 
b) various M and fixed n = 1 and = 0, c) two choices of n = 1, 5 and fixed M = 350 MeV and 
= 0. d) Axial form factor for different choices of M (there is almost no n dependence). 

VI. SUMMARY 

In this work we calculated analytically a set of photon distribution amplitudes up to 
twist four in tensor, vector and axial vector channels. We used nonlocal chiral quark model 
with momentum dependent quark mass. In order to get a correct behavior of low energy 
matrix elements we modified vector vertices (making them nonlocal) in such a way that 
Ward-Takahashi identities were fulfilled Similar, although numerical calculation was 
already done in ref. [l^. They also used instanton motivated nonlocal model with dressed 
vertices, taking into account rescattering in the p meson channel. The shape of the mass 
dependence on momentum was chosen as an exponent decreasing with k"^. Here we use F (k) 
as given by (flTl) and neglect rescattering which turns out to be small. 

First we obtained numerical estimates for quark condensate {i'lp), magnetic susceptibility 
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Xm. and decay constant f^^ in our model. For larger values of constituent quark mass M or 
power n our results are getting close to the ones of ref. {l^. Unlike (ipip) and Xm the value 



of /s^ is rather stable as far as model parameters are concerned. Using evolution equations 
(following refs. 1^) we find that for M = 350 MeV and n = 5, the scale of our model is 
about fi ~ 500 MeV (this estimation was done using {ipip), Xm, fs-y as given by sum rules 
at 1 GeV scale and evolving them backwards down to the model values). This is in rough 
agreement with the scale of the instanton liquid model which is believed to be of the order 
of 600 MeV 01 . 

Next let us discuss the properties of the distribution amplitudes obtained within the 
present approach. Leading twist amplitudes are not very sensitive to the value of power n. 
However it seems that higher twist DAs are rather strongly model dependent. 



Comparing our results with those of ref. [I8i| we find some similarities, but also some 
discrepancies. Tensor leading twist DAs are in fact the same. For real photons they are 
almost constant with small maximum at m = 1/2 and they do not vanish at the end points. 
The contribution of the nonlocal part of the vertex is rather small, it is however producing 
the small maximum in the middle. For ^ the end points move up for negative and 
down for positive P^, whereas the middle value behaves in the opposite way. 

Twist 2 vector DA (0^) should decouple for P^ = 0. Here the importance of gauge 
invariance shows up. We find cancelation of two contributions to 0y coming from the local 
and nonlocal parts of the photon vertex which are not proportional to P^. The remaining 
part is therefore proportional to P^ and decouples as required by the gauge invariance. We 
find that (pv vanishes at the end points and develops minimum in the middle for ^ 0, 
whereas for ^ it has a bell-like shape with a small dip in the middle. This behavior is 
very different from the one obtained in ref. {l^ where (py is almost fiat and does not vanish 
at the end points. However both vector and also tensor form factors are quite similar in 
both cases. 

One has to note that because of the subtraction of the perturbative part that is required 
in this case, (pv develops imaginary part for positive photon virtualities, so in this case we 
only discuss the real part. 

As far as higher twist DAs are concerned the situation is as follows. 

Our tensor twist 3 DA (i/jt) is identically zero for P^ = 0, because it is simply proportional 
to P^. For P^ < it is negative and has the shape of inverted "U", similarly to the one of 
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ref. 



18| . In this case DA is a regular function without S-tjpe singularities. Vector twist 



3 DA (ipv) in our case blows up at the end points; such behavior is not seen in ref, 
However, similarly to ref. j3| we also obtained delta-type singularities at the edges of the 
physical support. Twist 3 axial DAs (i/ja) in both cases show similar behavior: they do 
not vanish at the end points and have a minimum for u = 1/2. Despite the fact that for 

= axial vector DAs both in our case and in the case of ref. look similar, the axial 
form factors behave differently for < 0. In our case Fa^P"^) vanishes at large negative 
momenta in contrary to the one of ref. [l8| that tends to unity in the same limit. 

Regular part (without delta-type singularities) of twist 4 tensor DA {Iit) is in our case 



positive and vanishes at the end points for = whereas in ref. 



18| it is negative and 



does not vanish at the end points. For space-like photon momentum P^ < we see some 
similarity in shape between our hx and —hx of ref. fl^. Vector twist 4 DA {hy) is in our 
case a result of large cancelation of the positive non-local piece and the negative local piece. 
Its properties are not discussed in detail in ref. jisl . 

The only phenomenologically accessible photon distribution amplitude is the leading twist 
tensor DA - 0^. It is almost flat and does not vanish at the end points. This behavior is 
seen in our model and in other models discussed in ref. [18| and also in refs. Flat DA 

is characteristic for the elementary point-like particle, however, it is violating factorization 
theorems of QCD that require the DAs to vanish for n = 0,1. Formal evolution of such 
an amplitude is questionable, not only because the Gegenbauer series is not convergent at 
the end points, but also, because potentially large contributions coming from the vicinity of 



u 



0, 1 are not summed by the ERBL evolution equations 



3l|. 
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Appendix A: IDENTITIES 



In this appendix we summarize some of the identities used in this paper that deal with 
the sums of factors and powers of rji. Some of them have been already introduced in 
ref. but the general proofs have not been given. For definiteness let us recall that r]i, 
i = 1, . . . , 4n + 1 are the solutions of the algebraic equation G (z) = z^'^^^ + z^^ — = 0. 
We denote 

/i = n " ^3^^ ■ 

For any set of 4n + 1 complex numbers (not necessarily satisfying G [z) = 0) we have 

fl foTN = 4n, 
i=i I for iV < 4n. 

To prove flAll) let's define the following function 



M 

[li=i [z-Vi) 

where M < 4n and integrate it over a circle with infinite radius. On one hand we can use 
residue technique to get the sum entering (lAip . on the other hand, direct integration over 
the large circle gives right hand side of (lAip . 
If, in addition, rji satisfies G (rji) = then 

4n+l 

E /^^f = (-1)'^ (A3) 

for 4n < P < 8n. This can be proven in the following way. Notice that for P = 4n equality 
flA3[) is satisfied due to (lAll) . Let us move to P = 4n + 1, that is we want to calculate 

4n+l 

^ = E (A4) 

i=l 

Adding to this equation the result for P = An we get 

4n+l 

E (^'"''' + = ^ + 1- 

1=1 

-,4n+l I An Jl 



Using the fact that G irn) = rif^^ + ?7f" - = and ([Al]) we have 

471+1 

X + 1 = ^ = 0. (A6) 
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and X = —1. This procedure can be repeated several times to prove (IA3I) until P = 8n. 
For any set of 4?t, + 1 complex numbers rji we have the following identity 

4n+l 4n+l 



Imposing in addition the constraint G (?7i) = we get Vi = The proof is similar 

to the one of flAlj) . however we integrate the following function 

g (z) = 2 . (A8) 

{z - r]i) {z - r]2) ■ ■ ■ {z - r]i) ...{z- ?74„+i) 

After several algebraic steps we arrive at (jATj). 
Appendix B: PION DECAY CONSTANT 

In our model Birse-Bowler formula (fTSll for pion decay constant reduces to the following 
form 

- ' E Z^/^-^^l (1 + (1 + 2n)) r/f +i + (1 + 4n (1 + 3n)) r/f 



47^2 



+ 2n (1 + 6n) + 4n2r72n-2N /^^^ ^ - In (1 + 7],)) + (Bl) 

^ ^ \Vi-Vj ^ + ViJ 

where e^j is for i = j and 1 otherwise, while 6ij is Kronecker delta. 

Appendix C: LIGHT-CONE INTEGRALS IN SCHWINGER REPRESENTATION 

Here we summarize formulae used to perform d^k loop integration in the presence of 
6{n ■ K — un-p). We will consider three cases when the numerator contains no k~ at all and 
one or two powers of k~ . We follow closely the method of ref. jl^. 

Consider loop integral (l26l) and apply to it ([281) : 



expressed in terms of scaled variables (l2ll) . Here M is the numerator to be specified later. 
Recall that 

z^ = [ti-pf -l + ie, Z2 = -l + ie. (C2) 
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We shall now make continuation to the Euclidean metric: 



K° = iK^ (C3) 

with 



^ —K^, K-p^—K-p, n ■ K ^ —n ■ K, (C4) 



where arrows denote D dimensional Euclidean vectors. Therefore 



J {271 f ^ ''{i^^ + l + 7]i){{f^-p)^+l + 7]j) 



We shall parametrize now 



1 



oo oo 



{N? + l + r]i){{K-p)^ + l+Tlj) 

and 



-oo 



2s 



In these new variables we have 

471+1 



■i\p^ iu—y) 



(C5) 



da y rf/3e-"('^ +^+^^)-^(('^-*^ +^+^^") (C6) 





oo 

5{n-K + up-^) = f ^e~^Mn-^+up+)_ ^^^^ 



It is convenient to introduce new variables: 

a + P = s, P = ys, a = (1 — y)s = —ys. (C8) 
Integration measure reads then 

oo oo oo 1 

jdajdp = jsdsjdy. (C9) 



Finally we will shift momentum 

K = k' + ( yp- i—fi] . (CIO) 



I sdse-'V-y'^^+y'^^+yyA f ^Ll^j^e-'"^". (Cii) 

J J (27r) 



Further calculations depend on the nature of A/". If Af can be expressed entirely in terms 
of Zi^2 then, in virtue of (l28l) . it is enough to replace pertinent powers of Zi^2 ~^ Vi^j ^ind 
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perform Gaussian integration over k! . Let us denote such an integral as Jo- Also integral 
over d\ is trivial. In the following we shall need also integrals with A and which read: 

Hence (for D = A — 2e) we get 

= (4^ J ^ ^'"^^ /.^f j dss^-^-'V-^'^^^-^^^''^^ ] . (C13) 

In order to perform the integral over ds we shall use 

00 n £ 

/"dss^-i-"e-^[-l = [■■■f-'Vie-n) ~ — 1^^- -6"^^ + ... (C14) 

J e{e~l)...{e-n) 


arriving at 

^° 167r2p+ V ; £ ' [1 - m + uVj + «w P^T 

= 'l6^ ( ^ ) ^ ^^-^^'"^ ( - - In [1 - «r,. + uv, + uup^] V (C15) 
If numerator A/" involves additionally one power of k'^ we have then 

J\f ^ Afir]i, T],) {wk)= -J\f{7]^, 7]j) {w ■ k) (C16) 

where w is a constant four- vector. Let us denote such an integral by J7i. Here the only 
difference from the previous case comes from the integration over k . Since 

w ■ K = td ■ k' + ^ ■ p — i— w ■ fi^ (C17) 

only the terms in parenthesis survive. After integrating over dX with the help of flC12l) and 
over dy (in the case of 6' we have to integrate by parts) we arrive, back in the Minkowski 
metric, at: 

^ _ . A Avre-n ^ 1 _ ^ 



i,j=l J J J 

[u {wP)+ ((ry. - r/,) + (1 - 2u)P' 

[w ■ n) A^ 



2P+ e - 1 



(l + r7,)5(M-l) -{l + r]i)6{u) \ (C18) 
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where = P^/A^. Note that ii w = n then w ■ n = and we get J'q of eq. (ICISP 
multiplied by uP~^ as it should be, since we could have used 5 {k ■ n — un ■ P) in the first 
place. Similarly liw = e_[_ we have J7i = which means that a single power of k_l integrates 
to zero. Note that due to Lorentz invariance after du integration the coefficient in front of 
w ■ n should vanish in accordance with (l23l) . 

Finally if the numerator contains k^k*^, let's call such an integral J2, we have 



(C19) 



Using (1C12I) and integrating over dy we get three different contributions to J2 depending on 
the tensor structure: 

'I67r2p+ U7rA2 j e2^^'^Ml 

1 [1 — ur]i + ur]j + uup'^] 

2 (F^) 



J, 



(0) 



77^21 



urji + -ur^j + uup 



[wv) + — {vP){wP)\, 



(C20) 



^ 167r2 2P+2 V47rA2 ^ e 



(w ■ P) (f ■ + (if ■ ra) (t> ■ P) 



4n+l 



{(1 + 7^,)<5(W-1) 



and finally 



M?7i + ur]j + uup'^Y 
— (1 — urji + urjj + uup^) + m (77^ — r^^) + (1 — 2^)^^ 



4n+l 



(C21) 



(2) _ . ^ A^ / e-^ y " 1 



{t], - r]j) + (1 - 2M)p2] [[1 + .] 5(u - 1) - [1 + V^] ^iu) 

+ [1 — urji + urjj + uup'^~\ up^ + (r^j — r/j) + (1 — 2u)p'^ 



+ 



(C22) 



In eq. (IC22P we encounter derivatives of 5 functions; it is here implicitly assumed that the 
coefficient N [\ — urji + urjj + uup^] ^ when multiplied by — 1) or 5'{u) is taken at the 
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corresponding value of u. Note that Lorentz invariance requires that 

1 

j du J^^'^^ = (C23) 



(modulo possible subtraction of the perturbative part). 

Finally let us remark that if we need an integral of k'j_ we may use the following trick in 
two dimensional transverse plane: 

(fki_ kl = '^J d^k± {e± ■ k±f (C24) 

if there is no other dependence on the transverse angle, as it indeed happens in our case. 
We can then evaluate the r.h.s of eq. (IC24I1 using the formulae from the present appendix. 



Appendix D: AXIAL FORM FACTOR 

In this appendix we show, as an example, simple calculation of the axial form factor. We 
start from the matrix element on the left hand side in eq. (fT2ll : 

(0|^(-An)7^5^(An)|7(P,^)> = -eiVce. / gi(2/c.n-P.„)A ^j^^^ 

J (27r) 



Calculating the trace and taking the derivative with respect to A as in the definition (fT2i) 
we obtain 

Using Lorentz invariance and some simple algebra we get 

_ 9p/v p P „ / ^""^ {2k-n-P-n){k-fi-pk-n) 
M - -2eN^e.P^n,e J ^ D (k) D {k - P) ^^^^ 

with p = P^/P+^. Comparing this with the right hand side of (fT2il we get the following 
expression 

Fa{P') = ^-^J, (D3) 

where JT" denotes the integral in (|D2D . However it can be easily shown that Lorentz invariance 
requires that 



d^k k ■ n k ■ n — p [k ■ 2 f d^k k^ 



{2nf D{k)D{k-P) 2-Dj (27TfD{k)D{k-P) 



(D4) 



36 



and 

jD 



(2nfD{k)D{k-P) 

Using this, J' reduces to 



/d k k-n — pk-n /t\ \ 

7^ D{k)D{k~P) = 



4 r _<Pk kl 

^~2^j D{k)D{k-Py 

This integral reads 

aD-2 a- /I \D/2 4n+l i 

J=^^^[-^) E/^MNf / du{l-uv. + uv,+uup'Y-'r{e). (D7) 
Expanding in e and subtracting the perturbative part we obtain the final expression 

/ au [1 — urji + urij + uup ) m ^-^ . 

Jo V / 

Integration over du can, in principle, be done analytically (taking into account remarks given 

in the main text), however here we just plot Fa (P^) in sect. IV B 31 
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